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We analyze the above-threshold behavior of a mirrorless 
parametric oscillator based on resonantly enhanced four wave 
mixing in a coherently driven dense atomic vapor. It is shown 
that, in the ideal limit, an arbitrary small flux of pump pho- 
tons is sufficient to reach the oscillator threshold. We demon- 
strate that due to the large group velocity delays associated 
with coherent media, an extremely narrow oscillator linewidth 
is possible, making a narrow-band source of non-classical ra- 
diation feasible. 



Stable and low-noise sources of coherent and non- 
classical radiation are of interest in many areas of laser 
physics and quantum optics. Such sources have a wide 
range of applications such as frequency standards, optical 
magnetometry, gravitational wave detection, and high- 
precision spectroscopy. 

The present theoretical work is motivated by recent 
experiments demonstrating a phase transition to mirror- 
less oscillation of counter-propagating Stokes and anti- 
Stokes fields in resonant, double- A Raman media |Q. In 
contrast to earlier studies involving instabilities in alkali 
vapors (2]-[|] , this oscillation could be achieved with pump 
fields of fiW power (nano Joule pulse energy) and is ac- 
companied by a dramatic narrowing of the beat signals 
between driving and generated fields. Oscillations of this 
kind are a clear manifestation of atomic coherence and 
interference effects, which have recently lead to many ex- 
citing developments in resonant nonlinear optics p-p|. 
In particular, the unusual efficiency of the present pro- 
cesses is expected to lead to a new regime of quantum 
nonlinear optics in which interactions at a level of few 
light quanta are feasible. Furthermore, the photon pairs 
generated can possess nearly ideal quantum correlations, 
resulting in almost complete squeezing of quantum fluc- 
tuations [[To). 

We here study theoretically the quantum dynamics of 
the mirrorless oscillator above threshold. We show that 
for an infinitely long lived atomic dark state an arbi- 
trary small stationary flux of pump photons is sufficient 
to maintain the oscillation. Wc furthermore analyze fre- 
quency locking and linewidth narrowing of the beat note 
between oscillation- and pump frequencies. In particu- 
lar, we show that the beat-note linewidth is given by an 
expression similar to the Schawlow-Townes formula for 
lasers where the cavity storage time is replaced by the 



group delay time r gr in the medium. Due to the large 
linear dispersion associated with electromagnetically in- 
duced transparency (EIT) in optically thick media, the 
group delay can be extremely large 1 1 1 - 1 3 1 leading to a 
very small beat-note linewidth. This effect is analogous 
to the line- narrowing in intracavity EIT [pr|Jl5| . Since 
only very small pump powers are needed to reach thresh- 
old, ac-Stark shifts and the associated systematic effects 
on the beat-note frequency can be made very small. The 
combination of line-narrowing and small pump-power re- 
quirements makes the mirrorless parametric oscillator an 
interesting novel source of stable and narrow-linewidth 
non-classical radiation. Possible applications include fre- 
quency standards, optical magnetometry and few-photon 
nonlinear optics. 

Consider the propagation of four nearly resonant plane 
waves, parallel or anti-parallel to the z axis, in a medium 
consisting of double- A atoms (see Fig. 1). These include 
two counter-propagating driving fields with equal fre- 
quencies Ud and (complex) Rabi- frequencies Ef and Ej,, 
and two generated fields (anti-Stokes and Stokes) with 
carrier frequencies v\ and obeying v\ + v-i = 2ud- The 
fields interact via the long-living coherence (decay rate 
7o) on the transition between the ground state sub-levels 
bi and 62 with frequency splitting ujq = — 0Jbi ■ 

Due to resonantly enhanced four-wave mixing the 
coherent pump fields Ef and Ef, generate counter- 
propagating anti-Stokes and Stokes fields (here described 
by the complex Rabi- frequencies E\ and E2). For a 
sufficiently large density-length product of the medium 
and for a certain pump field intensity, the system shows 
a phase-transition to self-oscillations Q. The feed- 
back mechanism required for oscillations is provided here 
by the gain medium itself: A spontaneously generated 
Stokes photon stimulates "downstream" a Raman pro- 
cess. As a result an anti-Stokes photon is generated with 
a fixed relative phase. This photon propagates in the op- 
posite direction and stimulates another scattering event 
"upstream" . If the phase matching condition is fulfilled, 
this causes a second Stokes emission in phase with the 
first one closing the feedback cycle. A crucial condition 
for the coherence of this feedback mechanism is a suffi- 
ciently long-lived Raman coherence. 

We now discuss the transition to self-oscillation and the 
classical and quantum dynamics of the oscillator above 
threshold in detail. To simplify the analysis we ignore 
inhomogeneous broadening and assume equal coupling 
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strength of all four fields as well as equal radiative decays 
out of the upper levels. Furthermore we assume that the 
forward driving field Et is in resonance with the b 2 — * 0,1 
transition, whereas the backward driving field Ej, has a 
detuning A 3> \Eb\ from the b\ — > a 2 transition. In 
this case linear losses of the fields due to single-photon 
absorption processes are minimal. 

In order to calculate the medium response we solve 
the single-atom density matrix equations in third or- 
der of the Stokes and anti-Stokes fields and assume 
|A| > 7, \Ef t b\ > 7o,<5, where 8 = Vd + — v\ is the 
two-photon detuning. In a frame rotating with the car- 
rier frequencies, the propagation of the classical fields 
can then be described by the following equations for the 
slowly varying complex Rabi-frequencies 
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In these equations we have kept ac-Stark induced phase 
terms only in lowest order of the generated fields, since 
we are interested in in the case \Ef : b\ 3> |-Ei,2|- « = 
(3/87r)iVA 2 7 a is the equal coupling constant of all fields 
with N being the atom density, A the average wavelength 
of the fields, and 7 a the common population decay rate 
out of the excited states. E\ = E ie K^k- K /A)z ^ Ef> _ 
Ej,e tKZ ^, with Afc = k 2 — k\ being the phase mismatch. 

We note an important feature of the nonlinear coupling 
in Eqs. (0) and (0): In contrast to usual -media, the 
lowest-order cross-coupling terms are proportional to the 
ratio of the pump fields rather than the product. 
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Thus for \Ef\ = \Eb\ these terms are independent of the 
pump-field amplitudes. We will see later on that this 
leads to a rather unusual threshold behavior. 

In the present system a transition to spontaneous os- 
cillations is possible (!]], if the phase matching condition 
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is fulfilled. For large values of k, Eq.(||) describes a 
pulling of the frequency differences between generated 
fields and driving fields towards the ac-Stark shifted fre- 
quency of the Raman transition 
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This equation shows a close analogy with intracavity 
EIT. i] = cn/2\Ef\ 2 is a frequency stabilisation factor 
|l4| . This factor also governs the group velocity of the 
eigenmodes of the system v gr — c/(l + rj) and can be 
rather large. For conditions close to the experiments of 
Ref. fjjl, a reduction factor of 77 ~ 5 x 10 6 was mea- 
sured |13|]. In the limit of large rj the beat-notes between 
generated and pump fields locks tightly to the Raman- 
transition frequency of the medium. 

We next consider the classical steady state solution of 
the propagation problem. In the ideal limit (70 = 0) 
Eqs.([jJ|4|) have four constants of motion: the total inten- 
sity of the generated and pump fields |-Ei| 2 + \E 2 \ 2 and 
\E f \ 2 + \E b \ 2 , as well as Re \E* f E* h E x E 2 \ which has a simi- 
lar structure to the cubic expression conserved in 3-wave 
mixing Jr|. There is also the somewhat unusual con- 
stant of motion, l^/pexpdSil 2 /^ 2 ) 0- If we take 
into account, however, that Eqs. only hold to third 
order in the generated fields, this constant is equivalent 
to \Etf + \E x \ 2 . With this Eqs.(g-|) can be solved analyt- 
ically, if the phase matching condition is approximately 
fulfilled. Assuming equal input intensities of the driving 
fields 1 25/ (0)| = \E b {L)\ at z = and z = L respectively 
[L being the cell length), and disregarding linear losses 
due to the finite lifetime of the ground-state coherence, 
one finds in second order of the generated fields 



\E x {z)\ = E sini?(z), 



\E f (z)\=[\E f (0)\ 2 -E 2 sin^(z)_ 
\E b {z)\ = \\E f (0)\ 2 -E 2 cos 2 tf(z) 



\E 2 (z)\= E cos d(z), (8) 
1/2 

(9) 
(10) 



1/2 



where &(z) = Kz/A(l-E 2 /2\E f {0)\ 2 ). For kL/A < 
7r/2, E = 0. For values of kL/A larger than the critical 
value 7r/2 there are two solutions 



E = V2\E f (0)\ 
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with E — being unstable. It should be noted that in 
contrast to degenerate 4- wave mixing in usual x^ -m edia 
, the threshold condition does not contain the ampli- 
tude of the pump fields. Thus under the ideal conditions 
assumed here, i.e. for an infinitely long lived dark state, 
an arbitrarily small stationary pump intensity is suffi- 
cient to reach the oscillation threshold. Fig. 2 shows the 
above-threshold behavior of E as a function of kL/A and 
the field amplitudes normalized to |25/(0)| inside the cell. 
If the system oscillates not too far above threshold, the 
depletion of the pump fields is small and we may as- 
sume in the following constant driving field amplitudes, 
\E f (z)\ = \Eb(z)\=E d . 

To calculate the linewidth of Stokes and anti-Stokes 
fields relative to the drive field above threshold, we as- 
sume that the generated fields can be represented as 
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a sum of the classical stationary solutions and a time- 
dependent fluctuation (Ei >2 (z, t) — Ei >2 (z) + 8Ex, 2 {z,t)). 
We utilize a standard linearized c-number Langevin ap- 
proach in which collective atomic variables and fields are 
described by time- and position- dependent stochastic 
differential equations with S— correlated Langevin forces 
p9| . The diffusion coefficients or noise correlations are 
derived using the fluctuation-dissipation theorem and 
generalized Einstein relations. We obtain for the Fourier- 
components of the Stokes and anti-Stokes fluctuations 
8E x>2 {u) = 1/VZtt J dtSE h2 (t)e- luJt 
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where small fluctuation frequencies [uj <C Ed) and a con- 
stant phase of the pump field have been assumed. The 
term —kuj/E"^ is the contribution from the medium dis- 
persion. Following the procedure of Ref. we find for 
the dominant noise correlation 



(fx(z,uj)f 2 (z',uj')) 
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S{z-z')S(uj + u'), (13) 



where Af is the number of atoms in the cell, and we have 
identified the quantization length defined in Jl9| with the 
length of the cell L. 

Solving the inhomogeneous boundary problem for 
the Fourier-components of the field fluctuations with 
SE*(0,oj) = and 8E 2 {L,uj) — 0, we eventually find 
for the phase of e.g. SEi 
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In phase-diffusion approximation, the linewidth Ai>i of 
E x relative to the pump field is given by 



hx{L,w)54>x{L,J)) = Avi 



5{uj + u>') 



Using Eqs.([l3|) and (|l4|) we arrive at 



(15) 
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where P ou t is the output power of the mode. 

Eq.([l6|) can be represented in a very instructive form, 
if the group-delay time r gr = L/c(l + rf) in the resonant 
medium is introduced. In the appropriate limit, rj 3> 1, 
and near threshold such that kL/A « ir/2 the linewidth 
can be written as 
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Eq. ( |l7| ) is formally identical to that of an ideal laser with 
the cavity decay time replaced by the group delay time. 



In usual 4-wave mixing, based on non-resonant Kerr- 
nonlinearities , the group velocity is essentially equal 
to the vacuum speed of light. In the present scheme, how- 
ever, it can be substantially reduced due to EIT. 

It is important to emphasize that the photon pairs gen- 
erated by the oscillation process near threshold are in 
quantum- mechanically correlated states. This results in 
a dramatic suppression of intrinsic quantum fluctuations 
in a quadrature of the combined mode jl(| . 

In the discussion above we have neglected the relax- 
ation rate of the ground state coherence 70. If this de- 
cay is taken into account, one finds the modified thresh- 
old condition: cos(£L) + (7 A)/(2£^)sin(£L) = 0, with 
£ = k^I/ A 2 — 7q/4£^. In particular, oscillation can be 
achieved only if E 2 d > 7o|A|/2. This can be translated 
into a condition for the photon flux i.e. the number 
of pump-photons traversing the cell per unit time. One 
finds that the threshold photon flux in each pump beam 
is equal to the number of atoms in the ensemble decaying 
out of the dark state per unit time: 



$th = fM 70, 



(18) 



where / is a numerical prefactor of order unity. Since 
by using buffer gases or coated cells very small values 
of 70 can be achieved, a threshold flux corresponding to 
only few photons in the cell is feasible, leading to an 
interesting new regime of nonlinear optics. 

Furthermore, the non-vanishing linear losses resulting 
from the decay of the ground-state coherence lead to an 
additional noise contribution to the linewidth 
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This result can easily be interpreted. P ou tT gr /fiv is equal 
to twice the number of Stokes or anti-Stokes photons in 
the cell. (Note that the total number of Stokes plus 
anti-Stokes photons is spatially homogeneous.) Like in 
a usual laser, photon correlations are maintained over a 
time equal to the number of photons multiplied by the 
time a single photon stays in the system |^(| . The latter 
time is here given by the group delay time (if 70 is suffi- 
ciently small). If the lifetime of the dark state becomes 
shorter than the group delay, the phase information car- 
ried by a photon is lost faster and is dominated by 
270. Thus the minimum linewidth is ultimately deter- 
mined by the lower-level coherence decay. 

Similar to the case discussed in Ref. [Q for the intra- 
cavity system, the present result for the frequency lock- 
ing, Eq.(0), and the linewidth, Eqs.(16 1^) are a con- 
sequence of the large atomic dispersion associated with 
two-photon resonances in phase coherent media. In the 
limit of long-lived ground-state coherences, the beat-note 
linewidth can be extremely narrow. At the same time 
the resonantly enhanced nonlinearity makes it possible 
to achieve oscillation with very low pump powers. 

In order to see, whether the small intrinsic linewidth 
can indeed be exploited, we now estimate the influence 
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of systematic effects on the beat-note frequency. The 
most serious limitations arise from the ac-Stark shifts as 
indicated by Eq.(|^). At large values of pump intensi- 
ties these shifts are large and hence fluctuations in pump 
powers and frequencies will result in associated broaden- 
ing of the oscillator linewidth. However, the resonantly 
enhanced nonlinearity makes oscillation possible already 
when E\ > 70A, i.e. when the near-resonant ac-Stark 
shift E%/A exceeds the ground-state coherence decay 70. 
Thus with stabilized pump frequencies and intensities, 
technical fluctuations of the beat-frequency due to ac- 
Stark shifts could be several orders of magnitude smaller 
than 7o. In the experiment of Ref. jjj, for instance, short- 
term linewidth values below 100 Hz have been measured 
even though the transient time broadening of a Raman 
transition was about 50 kHz. It is clear that observa- 
tion of quantum-limited linewidth of the oscillator is most 
likely in the regime of ultra-low pump intensities. It is 
however this regime which is of main interest for quan- 
tum control and manipulation of quantum properties of 
few photon fields . 

In conclusion, we have demonstrated that resonant 
nonlinear interactions involving atomic coherence can be 
used for efficient generation of non-classical photon fields 
with a stable and narrow beat-note linewidth and small 
pump requirements. We expect these features to be of 
interest in many areas of quantum and nonlinear optics. 
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FIG. 1. Atoms in double A configuration interacting with 
two classical driving fields in forward (Ef) and backward di- 
rection (Et) and two quantum fields (£1,2). All optical tran- 
sitions are assumed to be radiatively broadened. 
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FIG. 2. Phase transition to mirrorless parametric oscilla- 
tions. Analytic solution for amplitude of generated field for 
70 = 0. Insert: Normalized field amplitudes inside medium 
for E/\E f (0)\ = 0.2 
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